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Abstract 

Cold stress causes various injuries to rice seedlings in low-temperature and high-altitude areas and is therefore an important 
factor affecting rice production in such areas. In this study, root conductivity (RC) was used as an indicator to map 
quantitative trait loci (QTLs) of cold tolerance in Oryza rufipogon Griff., Dongxiang wild rice (DX), at its two-leaf stage. The 
correlation coefficients between RC and the plant survival rate (PSR) at the seedling and maturity stages were -0.85 and -0.9 
(P = 0.01), respectively, indicating that RC is a reliable index for evaluating cold tolerance of rice. A preliminary mapping 
group was constructed from 151 BCjF, plants using DX as a cold-tolerant donor and the indica variety Nanjing 1 1 (NJ) as a 
recurrent parent. A total of 113 codominant simple-sequence repeat (SSR) markers were developed, with a parental 
polymorphism of 17.3%. Two cold-tolerant QTLs, named qRClO-1 and qRClO-2 were detected on chromosome 10 by 
composite interval mapping. qRClO-1 (L0D = 3.1, RMUI-RMUOS) was mapped at 148.3 cM, and qRCW-2 (L0D = 6.1, 
RM25570-RM304) was mapped at 163.3 cM, which accounted for 9.4% and 32.1% of phenotypic variances, respectively. To 
fine map the major locus qRClO-2, NJ was crossed with a BC4F2 plant (LI 88-3), which only carried the QTL qRClO-2, to 
construct a large BC5F2 fine-mapping population with 13,324 progenies. Forty-five molecular markers were designed to 
evenly cover qRCW-2, and 10 markers showed polymorphisms between DX and NJ. As a result, qRClO-2 was delimited to a 
48.5-kb region between markers qc45 and qc48. In this region, Osl0g0489500 and Osl0g0490100 exhibited different 
expression patterns between DX and NJ. Our results provide a basis for identifying the gene(s) underlying qRClO-2, and the 
markers developed here may be used to improve low-temperature tolerance of rice seedling and maturity stages via 
marker-assisted selection (MAS). 

Key Message: W\th root electrical conductivity used as a cold-tolerance index, the quantitative trait locus qRClO-2 was fine 
mapped to a 48.5-kb candidate region, and Osl0g0489500 and Osl0g0490100 were identified as differently expressed genes 
for qRaO-2. 
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introduction 

Rice [Oryza sativa L.) Ls one of the three major food crops in tlie 
world and feeds about 2 1 % of the global population in terms of 
per capita energy [1]. Cold stress is a common problem for rice 
cultivation and is a crucial factor affecting global food production 
[2] . A primary objective for rice breeding in low-temperature and 
high-altitude areas is to improve cold tolerance of rice cultivars. 
Therefore, mining cold-tolerant germplasms and mapping related 
genes would provide an opportunity to improve rice yield through 
marker-assisted selection (MAS) [3-5]. Recently, several quanti- 
tative trait loci (QTLs, e.g., Cth-1, Cth-2, qCTB3, qCTB7, and 
qCTBff) related to cold-induced rice damages to aboveground 



organs-including seedling mortality and growth, leaf wilting, and 
necrotic yellowing-were mapped [3,6-1 1]. To date, however, only 
Ctb-1, which improves the seed setting rate by enhancing spikelet 
fertility, has been isolated and cloned and has been proven to be 
involved primarily in the ubiquitin-proteasome pathway [8]. 
QTLs contributing to cold tolerance of seedlings are located on 
several rice chromosomes, accounting for 4-^9.3% of phenotypic 
variations [4,12-22]. In addition, certain chromosomes also bear 
germination cold-tolerant QTLs [23,24]. Cold tolerance of rice is 
controlled by multiple genes, and different cold-tolerance-related 
mechanisms take effect during the growth period of rice. 

Nevertheless, cold-tolerant QTLs have not been reported in rice 
roots, an important underground organ that protects rice from 
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cold stress. At low temperatures, the integrity of plant cellular 
membranes can be compromised by the accumulation of peroxide 
and oxide radicals, which increasing membrane permeability [25- 
27]. Rice roots are more sensitive to low temperature than leaves, 
and cold resistance at the seedling stage can be directly determined 
by the contribution of roots to cold resistance [28] . Therefore, root 
conductivity (RC) tests have been used to determine the degree of 
cold-induced injuries in wheat [29]. However, cold-tolerant QTLs 
have not been reported in rice roots based on the RC test. 

Dongxiang wild rice (DX), a Chinese common wild rice {Orpza 
nifipogon Griff, AA genome) grown injiangxi Province (1 16°36' E, 
28°14' N), is the northernmost wild rice known in the world [30]. 
It possesses extremely high tolerance to cold stress and is able to 
safely survive winter temperatures as low as -12.8°C [31]. In our 
present study, we reported fine mapping of a major QTL, qRClO- 
2, related to root cold tolerance using DX as a cold-tolerant donor 
and the cold-sensitive indica variety Nanjing 1 1 (NJ) as a recurrent 
parent. We also evaluated the effects of qRCl 0-2 at the seedling 
and maturity stages. The analysis of the QTL derived from the 
wild rice will provide novel alleles for improving rice cold tolerance 
by MAS and also genetic materials for understanding the 
molecular mechanism of cold tolerance in rice. 



Materials and Methods 

Plant Materials 

The cold-tolerant variety DX was used as a donor parent to 
cross with the cold-sensitive NJ in the summer of 2006 (Fig. 1). 
Then, a primary mapping population was constructed by 
employing the advanced backcross approach described by 
Thomson et al. [32]. The Fj progeny was first backcrossed to 
NJ to produce 151 BCiFi plants, which were then grown in a 
greenhouse, and subsequently each BCiFj plant was individually 
backcrossed to the recurrent parent until BC2F1 was obtained. No 
cold phenotype was selected for RC during population develop- 
ment because phenotype evaluation is destructive and thus cannot 
be conducted on individual plants. All BC2F2 seeds from 151 
BC2F1 plants were individually harvested, and the resulting 
population of 151 individual BC2F2 families was used for RC 
phenotypic characterization in September of 2008. 
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151 BC,F, plants X NJ 
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Figure 1. Process of constructing the mapping population. 

doi:10.1371/journal.pone.0096046.g001 
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Soaking time (min) 

Figure 2. RC of DX and NJ as the function of the time exposed 
to cold and soalting. A. RC measured as a function of cold treatment 
duration. B. RC measured as a function of soaking time. An asterisk (*) 
indicates that the RC between NJ and DX was significantly different (P< 
0.01). 

doi:10.1371/journal.pone.0096046.g002 

To eliminate the effect of qRClO-1, a single BC2F2 plant that 
contained the DX-introgressed segment at the locus qRClO-2 was 
selected for continuous backcrossing with NJ followed by MAS 
using the markers RM25570 and RM304. This yielded 14 BC4F2 
backcross inbred lines (BILs). Later, 150 simple-sequence repeat 
(SSR) markers covering the whole genome were used to determine 
the genetic background recovery rates, which varied from 89.33% 
to 94.67% (Fig. SI). Among these BILs, L188-3 (hereafter, L188) 
had a background recovery rate of 94.67%, which was found to 
have low RC. A large BC5F2 fine-mapping population of 13,324 
progenies was constructed by crossing NJ with LI 88 in the 
summer of 2012. Then, molecular markers were used to detect 
recombinant events covering qRClO-2. All needed recombinants 
were separately harvested and evaluated for RC, and the BC5F4 
seeds from the BC5F3 lines were harvested to calculate the plant 
survival rate (PSR) at the seedling and maturity stages. All plants 
were grown in the Yangzhou Wanton experimental fields at the 
Lixiahe Agricultural Research Institute of Jiangsu Province 
(119°42' E, 32°39' N) in summer and in Sanya of Hainan 
Province (110°02' E, 18°48' N) in winter. 

Evaluation of Cold Tolerance 

The RC phenotypes of the BC2F1 and BC5F2 plants were 
measured as the average RC of the derived BC2F2 and BC5F3 
plants across three replicates, respectively. Twenty seeds were 
randomly selected and placed in an illuminated incubator for 
growth to day 9. The seedlings were then given simulated diurnal 
alternating illumination treatment: 12 h dark at 4±1°C versus 
12 h with 25,000 Lx of illumination at 4±rC with relative 
humidity of 75-85%. After 48 h treatment, the 2-cm taproot 
section was immediately cut from each plant, and placed into an 
Eppendorf centrifuge tube with 2 mL deionized H2O (conductiv- 
ity of 1.5 |a.s/cm). After soaking for 1 h, the RC of each plant was 




2 3 4 5 6 7 8 9 10 ir 12 13 
RC (|is/cm) 



Figure 3. Cold-induced injury phenotype of parents and 
frequency distribution of BC2F1 plants after cold treatment. 

A. Cold-induced injury phenotype of DX and NJ after exposure to 4 'C 
for 48 h. B. Frequency distribution of the RC values among BC2F, 
population. Arrows indicate the mean scores of DX and NJ. 
doi:10.1371/journal.pone.0096046.g003 

measured with a conductivity meter (Mettler Toledo 326) and data 
were adjusted for RC of deionized H2O. 

Approximately 100 BC5F4 seeds were selected from each BC5F3 
line with high or low RC to evaluate the PSR using the criteria of 
RC<10 and RC>10 for cold tolerance and cold sensitivity, 
respectively. At day 9, the seedlings were cold-treated as described 
above and then returned to constant 28°C for recovery under 
illumination at 25,000 Lx for 72 h. The plants were divided into a 
sensitive group (with completely-withered leaves) and a tolerant 
group (normal growth). PSR was estimated as the proportion of 
the plants showing a tolerant phenotype to all the plants by 
percentage (0-100%). The PSR of each BC5F.5 line was measured 
in three replicates. 

Root cold tolerance at maturity was evaluated as follows: The 
BC5F4 seedlings obtained from the BC5F3 lines were planted in the 
field on May 10, 2013. Uniform seedhngs were then transplanted 
into a hydroponic greenhouse on June 10 at Yangzhou University. 
The hydroponic pool was 880 cm long, 130 cm wide, and 50 cm 
deep. The cement anchoring the rice plants was 135 cm long, 
16.7 cm wide, and 2.5 cm thick, and each had a plate of 14 holes 
that were 4.0 cm in diameter and away from each other by 10 cm. 
A single seedling was planted in each hole, which was fixed by cork 
and sponges. Each group (low RC vs. high RC) contained 126 
plants in three replicates. The hydroponic solution (pH 5.8) 
contained 10 p.p.m. of nitrogen (NH4NO3), phosphorus 
(NaH2P04), potassium (K2SO4), calcium (CaCl2), and magnesium 
(MgS04) each, 2 p.p.m. of iron (Fe(III)-EDTA), and 0.5 p.p.m. of 
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Figure 4. Genetic and pKiysical map covering qRCW-2. A. Location of qRC10-2 on rice chromosome 10. B. The linl<age map shows six indels and 
a SSR marl<er. C. Progeny test of homozygous recombinants delimited qRClO-2 to the region between the marl<ers RM25570 and RM304. The number 
of recombinants in each group and the phenotypic difference of each group from controls are shown. D. Further fine mapping of qRClO-2. E. Six 
genes in the candidate region. "Significantly different from PI at P = 0.01; ''significantly different from P2 at P = 0.01. PI: DX; P2: NJ. 
doi:1 0.1 371 /journal.pone.0096046.g004 



manganese (MnClj), and the solution was changed every other 
day. An oxygen pump was used to provide a steady flow of 
hydroponic nutrients, which maintained the nutrient concentra- 
tions and pH conditions in the pool during plant growth from 
seedling to maturity. These concentrations were elevated to twice 
the starting concentrations 60 days later. After 15-day post- 
heading, the cold hydroponic solution at 5±1°C was circulated 
around the roots for 24 h. After recovery at 28°C for two days, the 
PSR was measured. 



Linkage Mapping and QTL Analysis 

A genetic linkage map was constructed with Mapmaker/Exp 
3.0 [33]. The segregation data of each marker were statistically 
analyzed to measure the segregation distortion from Mendelian 
expectations. Table SI lists the SSR markers. QTL analysis was 
conducted by composite interval mapping using QTL Cartogra- 
pher 2.5 [34-36]. The logarithm of odds (LOD) threshold of 2.5 
was used to identify QTLs while avoiding false positives. 



Table 2. QTLs for RC detected in the BCjF, BlLs at seedling stage. 





Trait 


Chromosome 


Marker Interval 


Position (cM) 


LOD 


PVE 


Additive effect 


qRClO-l 


10 


RM171-RM1108 


148.3 


3.1 


9.4 


2 


qRClO-2 


10 


RM25570-RM304 


163.3 


6.1 


32.1 


4.3 



LOD: logarithm of odds; PVE: percentage of total phenotypic variance explained by the QTL; 
The LOD threshold for detecting QTLs by composite Interval mapping was 2.5. 
doi:l 0.1 371/journal.pone.0096046.t002 
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.2 s 
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Figure 5. Cold-induced phenotype of tKie BC5F3 lines contain- 
ing or lacking qRCJO-2 after heading at the 15-day stage. 

doi:10.1371/journal.pone.0096046.g005 
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DNA Extraction and Molecular Marker Analysis 

DNA was extracted from leaves following the cetyltrimethy- 
lammonium bromide (CTAB) method [37]. All leaves were stored 
at -80°C before extraction. Additional SSR markers between 
RM25570 and RM304 were chosen from the Gramene database 
(http://www.gramene.org). The sequences oi indica cv. 93-11 and 
japonica cv. Nipponbare were downloaded from IRGSPl.O (http:// 
rapdb.dna.afiFrc.go.jp/) and used to design insertion/deletion 
(indel) markers and real-time polymerase chain reaction (PGR) 
markers using Premier 5.0 (Tables S2 and S3). 

PGR and electrophoresis were performed following the method 
of Panaud et al. [38] . PGR was performed in a 20 |J,L volume 
containing 20 ng of template DNA, 0.15 |iL of 10 mM dNTPs, 2 
U of Tag DNA polymerase, 2 of 1 0 xPGR buffer (50 mM KGl, 
10 mM Tris-HGl pH 8.3, 1.5 mM MgGlj, and 0.01% gelatin), 
and 1.5 |J,L of 2 |J,M forward and reverse primers each. The 
cycling conditions were 5 min at 94°G, followed by 35 cycles of 
94°G for 1 min, 55°G for 1 min, and 72°G for 1 min, with a final 
extension at 72°G for 10 min. The PGR products were subjected 
to electrophoresis on 6% denaturing polyacrylamide gels, which 
were then stained with silver [39]. 



Data Analysis 

DifiFerences between parents and BILs for cold-related mea- 
surement were compared via Z-test using SPSS 1 1 .0 for Windows 
(SPSS Inc., 2002). The correlation coeflFicient between RG and the 
PSR was calculated using the SAS software (SAS Institute, 2000). 
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Gene Expression Analysis 

At the two-leaf stage after cold treatment for 0, 6, 12, or 24 h, 
the total RNA was extracted from 100 mg roots and 100 mg 
leaves using the RNeasy extraction kit (Invitrogen, Garlsbad, GA, 
USA), and any remaining genomic DNA was removed by DMasel 
treatment. Total RNA (3 |a,g) was used as a template for cDNA 
synthesis using M-MLV transcriptase (TaKaRa Biotechnology, 
Dalian, Ghina) with oligo (dT),^ primer. On an ABI PRISM 7500 
system (Applied Biosystems, Foster Gity, GA, USA), real-time 
PGR was performed in a 20 |a.L volume containing 2 |J.L of first- 
strand cDNA, 10 |XL of 2xSYBR Premix Ex Taq (TaKaRa, 
Shiga, Japan), 0.4 |J.L of 50xROX reference dye, and 2 |J.L of 
each primer (2.5 nM). The used amplification protocol was 95°G 
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Table 4. Candidate genes In the qRC10-2 region. 



ORF 


Putative protein function 


Pliysical location 


AA 


Gene size (bp) 


BAC ID 


Os10g0489448 


Unknown 


18,572,698-18,575,727 


125 


1425 


AC146481.2 


Osl0g0489500 


Terpene synthase 


18,572,698-18,575,857 


512 


1544 


AC 14648 1.2 


Osl0g0489650 


Unknown 


18,579,905-18,580,739 


101 


835 


AC146481.2 


Os10g0489800 


Unknown 


18,588,277-18,590,938 


384 


1363 


AC146481.2 


Osl0g0489900 


Putative DNA replication initiation protein 


18,591,175-18,595,632 


574 


2524 


AC 14648 1.2 


Osl0g0490100 


Pectin lyase fold family protein 


18,598,558-18,599,387 


214 


830 


AC146481.2 



ORF: open reading frame; BAC: bacterial artificial chromosome. 
doi:l 0.1 371 /journal.pone.0096046.t004 



for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 30 s. 
Rice Actin 1 was used as an internal control, and each experiment 
was repeated at least three times. 

Results 

RC is a Reliable Indicator of Cold Resistance in Rice 
Seedlings 

RC was significantly different for DX and NJ at day 8 (/■= 0.01, 
Table 1). The stem root was quite long and had excessive root 
hairs after 9 days, making it difficult to measure RC in an 
Eppendorf tube. Therefore the cold-treatment period was 
determined to be 9 days. The difference in RC between cold- 
sensitive and cold-tolerant plants was apparent for 1 h cold 
treatment (Fig. 2A). The RC for cold-tolerant seedlings increased 
to 8.1 at 24 h, but at 48 h the RC for DX decreased to 4.4 
whereas that for NJ increased to 12.9. In addition, NJ had high 
RC of 16.2 at 72 h whereas the cold-tolerant DX had low RC of 
5.3. The high RC in NJ can be explained by the observation that 
more intracellular liquids leak into the environment per unit time 
(Fig. 2B). 

QTLs for Cold Tolerance 

After recovery at 28°C, DX grew normally, whereas NJ showed 
cold-induced injuries (Fig. 3A). The RC was 2.9 for the cold- 
tolerant DX and 12.1 for NJ, and the distribution of the RC 
phenotypes among the BCjFi population was continuous and 
skewed towards tolerance (Fig. 3B). 

A total of 653 SSR markers were used to screen polymorphisms 
between the two parents. Only 113 markers were codominant, 
with the polymorphism of only 17.3%; all the 113 developed 
markers were mapped essentially and uniformly to the 12 
chromosomes. Based on this, two seedling cold-tolerant QTLs, 
named qRClO-1 and qRClO-2 (Fig. 4A), were detected on 
chromosome 10 via composite interval mapping. qRClO-1 
(LOD = 3.1) was mapped at 148.3 cM between RM171 and 
RM1108, and qRClO-2 (LOD = 6.1) was mapped at 163.3 cM 
between RM25570 and RM304, which accounted for 9.4% and 
32.1% of phenotypic variances, respectively (Table 2). 

Fine Genetic and Physical Mapping of qRCW-2 

To obtain a fine map of qRCl 0-2, 45 molecular markers were 
developed between RM25570 and RM304. As shown in Table S3, 
1 0 of these markers showed codominant polymorphisms between 
DX and NJ, consisting of 2 SSR markers {RAI25681 and 
RM25628} and 8 indel markers {qc8, qcl8, qcl9, qc25, qc27, qc45, 
qc48, and qc55). The arrangement of these markers is shown in 
Fig. 4B and 4D. 
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Figure 6. Expression of the candidate genes for plants 
subjected to 4°C cold treatment for 0, 6, 12, or 24 h as 
assessed by real-time PCR. 

doi:10.1371/journal.pone.0096046.g006 
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Using the markers RM25570 and RM304, a total of 21 1 BC5F2 
BIL recombinants were identified. Among them, 64 showed the 
homozygous genotype for DX or NJ alleles and were divided into 

six groups (LRl, LR2, HRl, HR2, Rl and R2) based on the 
genotypes of qc8, qcl8, qcl9, qc25, qc27, RM25681 and qc55 
(Fig. 4C). Group LRl (RC = 4.1) had DX alleles to the left otqc55, 
and group LR2 (RC = 3.5) had DX alleles to the right of 
RM25681, indicating that qRClO-2 was located upstream of qc55. 
Groups Rl and R2 had recombinant loci between RM25681 and 
qc55, and there was a large range in RC for these groups. Based on 
this, we deduced that qRClO-2 was located between RM25681 and 
qc55 (Fig. 4C). 

Based on the genotypes of the three polymorphic markers qc45, 
qc48 and RM25628 (Fig. 4D), groups Rl and R2 were further 
divided into four subgroups: Rl-1, R2-1, Rl-2 and R2-2. Then, 
according to the RC data, the candidate region was narrowed to 
68.7 kb between qc45 and RM25628, and qc45, qc48 and 
RM25628 were used to rescreen all mapped plants. Finally, we 
obtained four BC5F2 plants with high RC of 1 1.4 (named group 
R3) heterozygous in the qc48 allele. Therefore, qRClO-2 between 
qc45 and qc48 was restricted to 48.5 kb based on the Nipponbare 
sequence, in which a bacterial artificial chromosome clone 
(AC 146481.2) contained six candidate genes, as shown in Fig. 4D 
and 4E. 

Effects of qRClO-2 on Cold Tolerance at the Seedling and 
Maturity Stages 

Subgroups Rl-1, R2-1, Rl-2 and R2-2 were used to determine 
PSRs. The PSRs for Rl-1 and R2-1 with high RC values were 25- 
44% in seedlings and 21-35% in mature plants, whereas those for 
the subgroups Rl-2 and R2-2 with low RC values were 52-60'X) in 
seedlings and 48-57% in mature plants. The correlation 
coeHicients between RC and the PSR at the seedling and maturity 
stages were -0.85 and -0.9 (P=0.01), respectively (Table 3), 
indicating that a lower RC value was correlated with a higher 
PSR. Lines with qRClO-2 grew normally after 24 h under 5± 1°C 
treatment, and lines lacking qRCl 0-2 showed necrosis or death at 
the maturity stage (Fig. 5). These results indicated that qRCl 0-2 
contributed to rice cold tolerance at both stages. 

Identification of Cold-inducible Genes in the Candidate 
Region 

0. rufipogon Grifif, common wild rice, is the ancestor of cultivated 
rice [40^2]. The coUinearity of cold-tolerant genes between 
different species and tlie conservation in the cold-responsive 
pathway has been reported [30,43,44]. Therefore, Nipponbare 
(japonica) was selected as the reference genome to screen candidate 
genes. Table 4 lists the predicted functions of the six candidate 
genes in IRGSPl.O (http://rapdb.dna.afrrc.go.jp/). 

To further identify genes related to cold tolerance, real-time 
PCR was used to assess gene expression in the leaf and root tissues 
after 0, 6, 12, and 24 h of cold treatment (Fig. 6). The results 
suggested that Osl0g0489448, OslOg0489650, Osl0g0489800, and 
Osl0g04899()0 shared a similar gene expression profile in DX, 
L188 and NJ. Osl()g0489500 was induced under cold stress and 
showed high levels of inducible expression in the leaves of N^ at 
6 h and 12 h, whereas no significant variation was observed in 
DX. This pattern was different from that in the roots, as it was 
highly expressed in the roots of DX at 6 h but not up-regulated in 
NJ until 24 h. Osl0g0490100 in the roots of NJ showed higher 
expression than in DX before and after cold treatment. Therefore, 
Osl0g0489500 and Osl0g0490100 were considered as the genes 
associated with cold tolerance in the qRCl 0-2 region. 



Discussion 

Cold stress is one of the major environmental concerns in rice 
cultivation, especially in low-temperature and high-altitude areas 
[45] . The rice subspecies indica and japonica showed cold injury 
symptoms below 18°C and 15°C, respectively [46]. Under low- 
temperature stress, the root cell plasma membrane might become 
damaged and less fluid, resulting in functionality loss, which in 
turn caused increased cell membrane permeability and intracel- 
lular material's extravasation into the external environment 
[47,48,49]. In our present study, RC was used as a cold phenotype 
indicator reflecting the integrity and permeability of root 
membrane at low temperatures. With prolonged treatment time, 
the RC values of DX and NJ increased gradually, indicating 
cumulative cold injuries in both parents. Besides, more serious 
damages accumulated in the NJ root membrane. During the 
recovery period, high-RC plants showed wilting leaves and 
retarded growth, but low-RC plants maintained green leaves. 
Such a cold-tolerant phenotype had a consistent effect on plants at 
the maturity stage with qRCl 0-2. Moreover, the correlation 
coefficients (-0.85 and -0.9) between RC and the PSR were 
significant in seedlings and mature plants (P= 0.01), indicating that 
RC could be used as a reliable index for evaluating cold tolerance 
in rice. 

In previous studies, most cold-tolerant QTLs at the seedhng 
stage were detected using recombinant inbred lines derived from 
indica /japonica cultivars [12,18,50,51]. However, QTL analysis on 
a wider range of variations from wild species could be equally 
important. Koseki et al. identified the cold-tolerant QTL qCtsslO 
(LOD = 3.0) from W1943 (0. rufipogon) on chromosome 10, which 
accounted for 9.8% of phenotypic variations [17]. In our study, 
two QTLs qRClO-1 and qRClO-2 were identified on the same 
chromosome in DX. A QTL for cold tolerance (qCTSS-10) was 
also identified in seedlings by Yang et al. [52]. qCTSS-10 for 7°C 
cold treatment was mapped by high-throughput sequencing in the 
F3 mapping population from a cross between a cold-tolerant 
japonica rice variety (Nipponbare) and a cold-sensitive indica variety. 
The interval of qCTSS-10 (15.48-21.06 Mb) overlapped tiiat of 
qRClO-1 and qRClO-2, indicating that these QTLs may be 
controlled by the same cold-tolerant gene. Liu et al. also identified 
die same QTL {qCSTlO, LOD = 2.4) with 12% of the phenotypic 
contribution from the cold-tolerant line IL112 [53]; all early 
seedlings of the line ILl 12 survived normally for 9 days at 4— 5°C. 
A molecular marker tightly linked to qCSTlO was 0SR33, which 
was the same marker as RM171 linked to qRCl 0-1 and qRCl 0-2 
according to the Gramene marker data (http://www.gramene. 
org/). As described above, the same cold-tolerant QTL was 
identified on chromosome 10 among different parents and cold- 
treatment conditions based on a different cold tolerance index, 
suggesting that qRCl 0-2 was a stable and reliable cold-tolerant 
QTL for the rice seedling stage. Further studies on the candidate 
gene(s) in this region will aid in elucidating the molecular 
mechanism underlying rice cold tolerance. 

In the candidate region, Osl0g0489500 and Osl0g0490100 were 
differentially expressed between cold-tolerant and -sensitive plants. 
Osl0g0489500 encoded a terpene synthase that was localized to 
the cytoplasm. In our study, Osl0g0489500 was up-regulated in 
response to cold stress in the leaves of NJ. In the roots, however, its 
expression in DX and LI 88 was ~ 18-fold higher than that in NJ, 
suggesting that the root may be a key organ for the Osl0g0489500 
activity. Yang et al. also identified Osl0g0489500 as a gene up- 
regulated by low temperatures in a cold-tolerant variety [52]. After 
analyzing differentially-expressed transcripts in indica (K354) and 
japonica (C418) under controlled conditions (25°C) and cold stress 
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(4°C) for 48 h, Zhang et ad. identified two terpene synthase genes 
{Os02g36140 and Os08g07100) that were associated with signal 
transduction cascades fi-om K354 during cold-stress responses in 
seedlings [44] . Another gene differentially-expressed between cold- 
tolerant and -sensitive plants was OslOg0490100, which encoded 
virulence factors that degraded the pectin components of the plant 
cell wall [54,55]. Pectin functioned as an adhesive in cell wall 
formation, and the pectin content directly influenced cell wall 
integrity, which was related to rice cold tolerance as reported 
previously [56-58]. Solecka et al. found that cold acclimation 
induced an increase in the pectin content in the cell wall and 
increased tolerance to freezing [59]. A mutation in QUASIMO- 
DOl, a gene encoding the glycosyl transferase required for pectin 
s^Tithesis, caused reduced pectin content, leading to less cell 
adhesion and faster water loss during cold stress in common rice 
plants compared with the wild type [58]. Using Zhonghua 10 
plants, Liu et al. overexpressed OsBURPlG, a polygalacturonase 
ip that degraded pectin, and showed that an increase in the 
polygalacturonase activity led to a decrease in the pectin content in 
the cell wall, resulting in more severely curled leaves in the 
transgenic plants [53]. When the transgenic seedlings were 
exposed to 4°C for 72 h and allowed to recover for 7 days under 
normal conditions, the PSR of Zhonghua 10 plants was 42%, and 
that of the transgenic plants OE3 and OE4 was only 28% and 4%, 
respectively [60]. The above studies suggested that pectin 
degradation in the cell wall affected cell wall integrity and 
decreased toleranc;e to cold stress. In our present study, 
Osl0g0490100, a pectin lyase fold family protein that participated 
in pectin degradation, was more highly expressed in NJ roots than 
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